Comment on "Length-dependent translation of messenger RNA by ribosomes 
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In the recent paper of Valleriani et al [Phys. Rev. E 83, 042903 (2011)], a simple model 
for describing the translation of messenger RNA (mRNA) by ribosomes is presented, and an 
expression of the translational ratio r, defined as the ratio of translation rate cjti of protein 
from mRNA to degradation rate ujp of protein, is obtained. The key point to get this ratio 
r is to get the translation rate wti- In the study of Valleriani et al, wti is assumed to be 
the mean value of measured translation rate, i.e. the mean value of ratio of the translation 
number of protein to the lifetime of mRNA. However, in experiments different methods 
might be used to get wti. Therefore, for the sake of future application of their model to more 
experimental data analysis, in this comment three methods to get the translation rate Wti, 
and consequently the translational ratio r, are provided. Based on one of the methods which 
might be employed in most of the experiments, we find that the translational ratio r decays 
exponentially with the length of mRNA in prokaryotic cells, and decays reciprocally with 
the length of mRNA in eukaryotic cells. This result is slight different from that obtained in 
Valleriani's study. 



In rec- paper fl, V.lleri.ni a, presented a ,i„>- 

pie model to describe the length-dependent translation 
properties of messager RNA (mRNA). In their model, 
the mRNA degradation process is assumed to be gov- 
erned by rate lo^, i.e. the probability density of lifetime 
t of an intact mRNA is 



of translation rate uit\. Experimentally, there might be 
three methods to get ujt\: 

(I): LOti is obtained as the mean value of the measured 
translation rate f{t) from mRNA to protein, i.e.. 



^ti = {f{t)) = / f{t)<l>u[t)dt, 
Jo 



(3) 



= Ur exp(— w,.f). 



(1) 



Meanwhile, the rate of ribosome entering the coding re- 
gion of mRNA is assumed to be Won, and the degradation 
rate of protein is denoted by cjp. 

To discuss the mRNA length-dependent properties of 
the translation to protein by ribosomes, in jlj the expres- 
sion of translational ratio, defined as 



where f{t) = N{t)/T{t), N{t) is the mean number of 
proteins that the mRNA will synthesize if degradation 
occurs at time t, and T{t) is the lifetime of a mRNA 
before completely degraded {t is the lifetime of an intact 
mRNA). 

(II): ujt\ is obtained as the ratio of mean number {N{t)) 
of proteins translated from one mRNA to the mean life- 
time {T{t)) of a mRNA, i.e., 



(2) 



(iV(i))/(T(i)). 



(4) 



is obtained for translations in both prokaryotic and eu- 
karyotic cells. Where Wti is the translation rate from 
mRNA to protein by ribosomes. Since the protein degra- 
dation rate ujp is independent of mRNA, the essential 
point to analyze the mRNA length-dependent proper- 
ties of the translational ratio r, is to get the expression 



(III) : wti is obtained as the reciprocal of the mean dura- 
tion time of translating one protein from mRNA, 
1 1 



(l//(i)) {T{t)/N{t))- 
One can easily show that, for the mRNA transla- 
tion problem discussed in T{t) = t, N{t) = 9{t - 



2 



i^™)a;on(^ — tY°) ^'^^ translation in prokaryotic cells, and 
T{t) = t + fj^, N{t) = Uont for translation in eukary- 
otic cells. Where 9{t) is Heaviside function, i.e. 9{t) = 1 
for i > and e{t) = for t < 0, = L/vP''° and 
— are the time taken by ribosonies to reach 

the end of niRNA, L is the length of mRNA, yP™ and 
w"" are the average velocities of ribosonie along niRNA 
in prokaryotic and eukaryotic cells respectively. 

Intuitively, method (I) is reasonable, and actually this 
method is used in [l| . On the other hand, method (III) is 
usually employed for some mathematical problems. For 
example, to get the mean translation rate of a process 
which includes two sub-processes with rate uji and W2, 
one usually does the following calculations 



0J1UJ2 



(6) 



(T) (Ti) + (T2) l/ui + 1/^2 t^l + ^2 

It should be pointed out that, if the method (III) is 
employed to get luh, the experimental samples with no 
protein synthesized, i.e. samples for N{t) = 0, should be 
discarded to avoid the infinite waiting time cases. Cor- 
respondingly, in the theoretical calculation in Eq. (5), 
the average {T{t)/N{t)) should be done for only large 
enough time t > tum which satisfies N{t) > 1, and 
the probability density 0^/ should be changed accord- 
ingly, (t)uit) = (t>u{t)/ J^^ (t>u{t)dt = LOr e-xp[-ujr{t-t\im)\ 
for t > iiini. In our numerical calculations, we use 
^\im ^ + l/i^on for translation in prokaryotic cells, 
and tJim = l/'^on for translation in eukaryotic cells. 

Meanwhile, the simple method (II) is often used in 
the experimental data statistics. For this mRNA trans- 
lation problem, one can easily get the following results by 
method (II). For mRNA translation in prokaryotic cells. 



=Wo„exp(-a;rt£ ). 



(7) 



So the translational ratio rP™ = uj^° juj^ decays expo- 
nentially with the length LP™ of mRNA. Meanwhile, for 
the mRNA translation in eukaryotic cells, 

= (^o„i)/ (t + tT) - C^on/ (1 + ^rtt)- (8) 

So, roughly speaking, the translational ratio r'^^ — 
uf^ jujp decays reciprocally with the length of mRNA. 
Note: similar as the above discussion about the calcula- 
tion of translation rate wti in prokaryotic cells by method 
III, if the samples with N(t) — are discarded in exper- 
imental measurements, the methods I and II should also 
be changed correspondingly, with a modified probability 
density 0(7 (i) for t > tn^ only. 

For the sake of comparison, the theoretical results of 
translational ratio r of E. coli and S. cerevisiae, obtained 
by the three methods, are plotted in Fig. [TJ with model 
parameters listed in Tab. HI Since method I is used by 
Valleriani et al. in [l| , the curves in Fig. [1] (and the cor- 
responding model parameters in Tab. HI are copied from 
1]. One can see that, compared with the experimental 
data, these three methods all work reasonably well, but 
with slightly different model parameter (see Tab. |l|. The 
most reasonable parameters should be the ones which are 
obtained by the same method as that used in the exper- 
imental measurements. 

In conclusion, when one tries to apply the model pre- 
sented in [l| to real experimental data, the method used 
to get the translation rate Wti should be chosen properly 
to be consist with experiments. 
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TABLE I: Model parameters LUonUop and Ur obtained by fitting the theoretical results to experimental data of E. coli and S. 
cerevisiae respectively (see Ref. [l| and references therein for the detailed description of the experimental data) . 





E. coli 


S. cerevisiae 






t^on / 


a;r(min ^) 


Valleriani et al. 


708.2 0.2 


7.69x10^ 


0.2 


method II 


734 0.61 


6.34x10^ 


0.22 


method III 


841.1 0.15 


6.62x10^ 


0.02 
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FIG. 1: (Color online) Theoretical results of the translational ratio r^"^" for E. coli (a) and for S. cerevisiae (b). The dots 
and solid curves are copied from Ref. [H (see the related references therein). The model parameters LOou/ujp and ujr used in the 
calculations are listed in Tab. [H 



